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ABSTRACT 

As  efforts  continue  to  effectively  utilize  the  microwave  frequencies 
a  need  is  apparent  for  a  traveling  wave  type  device  capable  of  handling 
greater  power  at  higher  frequencies .     One  possible  solution  is  the 
plasma  microwave  amplifier  tube,, 

This  paper  discusses  an  analysis  of  the  plasma  -  helix  coupling 
problem  which  may  hopefully  lead  to  development  of  such  a  device 
whereby  the  R.  F.  energy  is  coupled  directly  to  the  plasma  with  a 
helix,  rather  than  coupling  to  an  electron  beam  as  has  been  the  case 
heretofore.    A  theoretical  analysis  is  presented,  and  the  computer 
solution  to  the  analytical  equations  evolved  and  the  resulting  oj-  j3 
plots  are  included. 

This  writer  wishes  to  express  his  appreciation  for  the  assistance 
and  inspiration  provided  by  Dr„  Glen  A„  Gray  of  the  U.   S,  Naval 
Postgraduate  School  in  carrying  out  this  investigation. 
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1.    Background. 

The  purpose  of  this  investigation  was  to  study  the  problem  of 
coupling  electromagnetic  energy  between  a  slow  wave  structure  and  a 
plasma  column. 

The  underlying  reason  for  interest  in  this  problem  is  the  design 
of  a  plasma  microwave  amplifier  tube.    A  plasma  microwave  amplifier 
tube,  where  the  plasma  column  acts  as  the  slow  wave  structure,  inter- 
acting with  the  electron  beam,  may  well  provide  some  advantages  over 
present  day  types  of  traveling  wave  tubes .    At  the  present  time,  con- 
ventional traveling  wave  tubes  are  limited  in  the  power  they  can  handle 
at  high  microwave  frequencies.    As  frequency  increases,  it  becomes 
difficult  to  construct  a  satisfactory  slow  wave  structure  which  will 
dissipate  any  appreciable  heat.    Even  at  the  lower  microwave  frequencies, 
heat  dissipation  in  the  slow  wave  structure  may  be  a  problem  for  high 
power  applications.    The  plasma  as  a  slow  wave  structure  offers 
definite  possibilities  in  that  heat  dissipation  in  the  structure  should  no 
longer  be  a  problem.     There  should  be  no  inherent  limitation  insofar  as 
mechanical  dimensions  are  concerned,  as  the  frequency  is  increased. 
Further,  the  plasma  microwave  amplifier  device  appears  to  offer 
promise  of  high  gain  and  high  power  in  a  relatively  small  package. 

The  device  proposed  (Fig.   (1) )  is  in  many  respects  similar  to  a 
conventional  traveling  wave  tube.    The  tube  consists  of  a  conventional 


-fH 


►  ■ 


electron  gun  assembly,  a  d.c.  heated  filament  consisting  of  a  .001" 
tungsten  ribbon  approximately  1/2"  X  3/4"  at  either  end  of  the  interaction 
region,  and  two  helix  slow  wave  structures  in  the  interaction  region. 

The  concept  is  to  generate  a  plasma  by  heating  the  two  filaments 
to  2300   C  in  a  cesium  gas  environment  at  a  very  low  pressure,  on  the 

-6  -7 

order  of  10       to  10       mm.  Hg.     These  hot  filaments  will  then  cause 
contact  ionization  of  the  cesium  gas  thereby  generating  a  plasma. 
R.  F.  energy  is  then  coupled  from  the  helix  slow  wave  structure  to  the 
plasma,  exciting  the  plasma.     In  turn,  the  electron  beam  on  passing 
thru  the  plasma  interacts  with  the  wave  excited  on  the  plasma,  giving 
a  traveling  wave  type  interaction. 

A  helix  slow  wave  structure  was  chosen  for  several  reasons.     First, 
it  is  rather  easily  fabricated.     Second,  it  has  been  used  extensively 
in  tube  work,  and  experience  is  available  regarding  its  characteristics, 
Finally,  a  theoretical  analysis  of  the  helix  is  available,  Beck  (1),  and 
H  utter  (2). 

The  cesium  plasma  was  chosen  for  several  reasons.    First,  the 
plasma  frequency  may  be  easily  varied  by  control  of  the  envelope 
temperature,  and  thereby  the  gas  pressure.     Second,   it  provides  a 
relatively  homogenous  plasma,  and  may  be  quite  highly  ionized.     Third, 
it  provides  a  quiet  plasma,  that  is,  it  is  relatively  free  from  spurious 
oscillations  and  noise.    Finally,  it  may  be  readily  ionized  by  contact 


ionization.     Since  the  ionization  voltage  of  the  cesium  is  below  the 
work  function  of  the  heated  tungsten  filament,  the   filament  will  emit 
electrons  and  at  the  same  time  ionize  the   cesium  gas,  generating  a 
plasma  which  is  relatively  homogenous  „     Knechtli  and  Wada  (7)  have 
done  considerable  work  on  this  type  of  plasma  generation.     Kino  (8) 
has  also  used  a  similar  scheme.     In  his  tube,  Fig.,    (3),  a  spiral  heater 
was  used,  and  due  to  the  fields  set  up  across  the  spiral,  it  was 
necessary  to  make  all  measurements  during  the  time  when  power  was 
not  applied  to  the  heater.    By  use  of  a  flat  ribbon,  it  is  hoped  to  be 
able  to  reduce  the  field  of  the  heater  sufficiently  so  that  measurements 
may  be  taken  continously  while  applying  ad.c.  voltage  to  the  heater. 

Work  on  plasma  microwave  amplifier  tubes  has  been  carried  on 
by  several  authors,   (8),   (9),   (12),   (13),   (14),  and  (15).     The  major 
difference  between  the  operation  of  the  tube  proposed  here  and  the 
work  of  some  previous  authors  is  in  the  manner  of  coupling  the  energy 
into  the  device.     Boyd,  Field,  and  Gould,   (9)  studied  plasma  oscillations 
and  growing  plasma  waves  using  a  mercury  arc  discharge  for  a  plasma, 
coupling  the  input  R.  F.  signal  to  the  electron  beam,  and  thence  to  the 
plasma.     Fig.   (2)  is  a  sketch  of  the  device  they  used,    Allen  and  Kino  (8) 
studied  the  interaction  of  an  electron  beam  with  a  fully  ionized  plasma 
using  a  cesium  plasma  and  coupling  the  input  R.  Fo   signal  to  the  electron 
beam,  and  thence  to  the  plasma.    Fig.   (3)  is  a  sketch  of  the  device  they 
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used.     Paik  (12)  used  a  coupled  mode  approach  to  study  the  coupling  of 
modes  between  a  slow  wave  mode  and  a  helix.    The  experimental  work 
he  reported  on  was  done  with  a  device  as  shown  in  Fig*    (4)  which  used 
a  mercury  plasma  „     He  excited  only  one  mode  by  careful  coupling, 
using  a  selective  mode  coupler.    This  device  did  not  include  an  electron 
beam,  and  his  theoretical  and  experimental  work  was  limited  to  the  case 
of  no  magnetic  field „ 

B.M.   Bulgakov,  et  aL    (13)  studied  slow  waves  in  a  helix  -  plasma 
column  with  an  arbitrary  magnetic  field  .    This  did  not  include  any  ex- 
perimental work,  and  the  final  dispersion  plots  they  developed  showed 
discontinuities  in  the  curves,  which  one  would  not  expect  to  see. 
This  has  caused  some  workers  in  the  field  to  question  their  results. 

E.  V.  Bogdanov,  eta. I.   (14)  studied  very  nearly  the  same  problem 
as  Boyd,  Field,  and  Gould  (9).    The  device  they  used  for  experimental 
work,  although  physically  different,  used  an  arc  discharge  type  of 
plasma  generator  and  coupled  the  input  R„  F,  energy  to  an  electron 
beam,  and  the  output  R.  F.  energy  from  the  electron  beam.    The  signifi- 
cant thing  about  this  work  is  they  were  able  to    obtain  gain  from  25  cm. 
up  to  3  cm.     This  was  possible  because  they  used  an  arc  discharge  to 
generate  the  plasma,  and  by  increasing  the  arc  current  could  generate 
a  dense  plasma,  which  was  then  contained  by  a  magnetic  field. 

Work  on  an  experimental  X  band  tube  has  been  earned  out  by 
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1 
Microwave  Associates,  Burlington,  Massachusetts,       This  is  a  mercury 

device,  again  coupling  the  signal  via  an  electron  beam.    Eicon  Labor- 
atory, Cambridge,  Massachusetts  is  working  on  a  device  using  the  fast 
wave,  designed  to  utilize  the  negative  conductance  of  the  plasma.    All 
other  devices  proposed  to  date  have  used  the  slow  wave  mechanism. 
General  Electric  Ltd.  ,  Wembly  England,  has  reported    a  plasma  X  band 
unit  which  develops  60  db  of  gain  in  a  six  cm.  tube.     Philips  Research 
Laboratories,  Eindhoven,  Netherlands  has  reported    a  helix  coupled 
tube  giving  19  db  gain  at  four  Gc. 

In  all  of  the  preceding  devices  coupling  of  the  input  signal  to  the 
device  was  accomplished  by  coupling  as  in  a  conventional  traveling 
wave  tube,  to  an  electron  beam  with  a  helix,  and  then  using  the  wave 
propagated  on  the  beam  to  couple  to  the  plasma.    As  a  result  of  this, 
in  the  experimental  work  performed  to  date  it  has  not  been  possible  to 
couple  directly  to  a  backward  wave  plasma  mode.     It  is  anticipated  that 
with  the  experimental  device  proposed,  Fig„   (1),  it  will  be  possible  to 
investigate  backward  wave  modes  by  direct  coupling,  as  shown  in 
Fig.   (5). 


M.  Meisels,  Plasma  Amplifiers  Move  Closer  to  Hardware  Stage 
Microwaves,  pp  4-6,  Jan.  1963 
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If  energy  is  coupled  into  and  out  of  the  plasma  in  this  fashion,  then  jf 
interaction  with  the  electron  beam  is  observed,  it  must  necessarily  be 
interaction  with  a  backward  wave.  Of  course,  care  must  be  exercised 
to  ascertain  whether  this  interaction  occurs  in  the  helix  or  in  the  drift 
region. 

In  considering  the  specific  problem  of  coupling  electromagnetic 
energy  to  a  plasma  column,  it  is  perhaps  worthwhile  to  review  the  work 
done  in  studying  propagation  on  an  electron  beam,    Hahn  (10)  and 
Ramo  (11)  considered  a  drifting,  ion-neutralized  cylindrical  electron 
beam  in  an  infinite  axial  magnetic  field  and  found  that  in  general  there 
are  two  space  charge  waves,  the  propagation  equation  being0 

0Jp2  -  (w  -  j8uo)2=  0  (1.1) 

for  the  one  dimensional  case,  where 
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e 

m 

and  all  quantities  are  assumed  to  have  e         "^      dependence .    Thus, 

the  solutions  are; 

CO  ±  cop 

p  =  — -  (1'3) 

If  Uq  is  zero,  this  one  dimensional  ion-neutraiized  electron  beam 
does  not  propagate  space  charge  disturbances.,     Trivelpiece  (6)  has 
shown  however,  that  a  finite  u0    is  not  essential  to  the  propagation  of 
space  charge  disturbances  when  the  beam  has  a  finite  transverse  geo- 
metry. 

Rigrod  and  Lewis  (3)  considered  the  problem  of  an  electron  beam 
under  Brillouin  flow  conditions  and  found  two  types  of  space  charge 
wave  propagation,  one  involving  a  perturbation  of  the  average  charge 
density,  and  the  other  surface  rippling.     Brewer  (4)  considered  the  more 
general  case  of  an  arbitrary  magnetic  field. 

2.    Statement  of  the  problem. 

In  order  to  study  the  coupling  of  R.F.  energy  between  a  helix  slow 
wave  structure  and  a  plasma  column,  it  is  desirable  to  determine  the 
to  -  j8  plot  or  Brillouin  diagram  for  the  proposed  circuit.    This  provides 


11 


all  the  information  required  to  forecast  the  performance  in  interaction 
with  electron  beams,  except  knowledge  of  the  actual  field  amplitudes 
involved,  ref.  Beck  (1)„ 

The  term  plasma  as  used  here  denotes  a  partially  ionized  gas 
which  is  electrically  neutral  in  any  macroscopic  volume  in  the  absence 
of  disturbances.,    Further,  for  this  discussion  it  has  been  assumed  the 
ions  are  stationary,  and  the  electrons  have  no  random  velocities  and 
suffer  no  collisions  „    The  neutral  gas  molecules  are  assumed  to  play 
no  part  in  the  process  studied „     It  is  felt  this  is  a  reasonable  model 
for  the  plasma.    Fig,   (6)  shows  the  basic  configuration  of  the  problem* 
The  helix  slow  wave  structure  has  been  appproximated  by  the  sheath 
helix  model. 

3.    Analytical  Solution . 

Trivelpiece  (6),  has  solved  the  problem  of  a  plasma  filled  cylindri- 
cal waveguide  in  an  infinite  axial  magnetic  field,  starting  with  the 
basic  Maxwell  equations,,     He  further  considers  the  finite  magnetic 
field  case  where  he  neglects  the  aoC  magnetic  fields,   (qua  si- static 
approximation) .     He  presents,  for  the  finite  magnetic  field  case,  the 
development  of  a  complex  dielectric  tensor,  to  take  into  account  the 
effects  of  the  plasma  in  a  finite  magnetic  field „    This  Is? 
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(3.4) 


Making  the  quasi-static  approximation,  i.e.  setting  the  a,c„  magnetic 
fields  to  zero,  we  may  solve  the  problem  of  matching  between  the  helix 
and  plasma  column  as  follows; 


7  x  E  =  -jccB  =  0 

:.  e  =  -v  0 

This  approximation  is  presumed  good  for  V    phase  <  <c 
Now,  for  the  equivalent  dielectric, 

7  •  D  =  7  °   (^    ■    E  )  =  0 
which  suggests  the  differential  equation 

7  .   _^     .     70   =  0 
The  solution  to  which  in  cylindrical  coordinates  is- 


(3.5) 


(3.6) 


(3.7) 
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In  general  one   may  assume  a  solution  of  the  form 

#-»W.-Me'*  (3.9) 

However,  we  shall  restrict  ourselves  to  the  case  of  no  variation  in  6 
(i.e.  n  =  0) .    This  is  not  essential,  but  will  simplify  the  analysis. 
Also  this  is  the  case  of  greatest  practical  importance.    This  then  leads 
to  the  partial  differential  equation 

r      dr  dr  €  d    2 

rr  z 


(3.11) 


and  solve 

1       d     .     dR    . 
r     dr     [T   dr    ' 

-^    €ZZ    R=0 
rr 

and  now  let 

rr 

^     1  ~  (WpVce2) 

*      l  -  [  oyV  Ceo2  - 

cc2)  ] 

(3.12) 


The  solution  to  this  equation,  then  is.  0  <  r  <  b 


E,  =  AJ0  (rr)e 


J<«€-/*z>  (3.13) 
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E    -  -A  --^-   Jx  |rr)e  *z 

r  P 


H0=-A—j2=-    Ji  (rr)e  *V 


(3.13) 


and  for    b  ^  r  ^  a 
E 


=  [BI0(iSr)+CK0(J8r)]e-{COt'  P  z] 


Er  =  j[Bl!(i8r)  -  CKX  (j3 r) ]  ej (0i t  "  ^z} 


(3.14) 


H0  = 


BUjSr)  -  CKx  (0r) 


licet-  j8    ) 
e  'z 


and  with  suitable  matching  at  the  boundaries,  we  may  reduce  the  above 
to° 


crr      Mr  a)    = 
0a       J0(Ta) 


(3.15) 


Ii(j3a)  [KoOSbl  +  FKxflSb^+frQSaKl, 


Fl^jSb)] 


-I0(j3a)[Ko(i3b)  +  FK!(/5b)}f  Kq  (/3a)[Io  (|8b)  -  FlxCjSb)] 

where  F  is  a  function  of  the  admittance  Y    of  the  sheath  helix,  defined 

s 

as; 


]8Y 


13.16) 


Y    is  the  admittance  of  the  helix.    The  expression  for  Y    is  given  by 
s  s 

Hutter  (2)  as: 
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-J  OLia    In(iSa)    .   KoOa)        Ki  (j3a)    ,ka_  3        - 

|8       h  (0a)      Kx  (]3a)        Ko  (j3a)  l0a  j    COt    * 
Y    = ____ (3    17) 

(^)2cot=* 

using  the  sheath  helix  approximation .     It  should  be  noted  that  inherent 
in  use  of  this  expression  for  the  admittance  of  the  sheath  helix  is  the 
further  assumption  of  a  free  space    TE    mode  within  the  heJix,    The  use 
of  a    TE    mode  in  addition  to  the    TM    mode  is  required  to  match  the 
boundary  conditions  at  the  helix. 

Paik,  (12)  has  developed  an  analysis  using  the  coupled  mode 
approach  for  a  helix-plasma  column  in  the  zero  magnetic  field  case. 
The  above  general  expression,  equation  (3  ,15)  may  be  shown  to 
reduce  to  that  of  Paik  in  the  zero  magnetic  field  case,   (let  co    go  to 
zero) . 

Fig.   14  shows  an  00-/3  plot  for  the  zero  and  infinite  magnetic 

field  case,  with  the  parameters  defined  as  follows: 

Ci     =    k  a  cot  ij) 
1  P 


c2 

=    b/a 

c3 

=    k  a  cot  lb 
c 

X 

=    pa 

y 

=    ka  cot  $ 

where  k  =  w/C 

C  = 

the  velocity  of  light 
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The  data  was  obtained  by  use  of  a  digital  computer.     The  general 
program  and  necessary  function  subroutines  are  given  in  appendix  I. 

Fig.   (7)  and  Fig.   (9)  thru  Fig.   (17)  show    cc  -  y3  plots  for  various 

parameters  in  the  general  magnetic  field  case.     Fig.   (7)  shows  the 

principle  modes,  as  well  as  a  few  of  the  higher  order  modes.     These 

higher  order  modes  result  mathematically  from  the  multiple  roots  of  the 

Bessel  functions,  and  represent  higher  order  modes  of  operation. 

Since  there  are  an  infinite  number  of  roots  to  the  Bessel  functions, 

there  are  an  infinite  set  of  these  higher  order  modes  in  the  regions 

where  the  propagation  constant  is  real.    A  sample  plot  of  this  situation 

is  shown  in  Fig.   (8)  for  a  typical  case.    The  symbols  used  follow  those 

of  the  computer  program.    Z  is  the  right  hand  side  of  Eq  „    (3  „  15)  and  BB 

is  the  left  hand  side  of  this  equation.    This  plot  is  for  a  fixed  value  of 

x.     In  the  region  where  the  propagation  constant  is  imaginary.  Z  reduces 

to  a  function  of  two  modified  Bessel  functions,  usually  represented  as 

I   (x)  which  do  not  exhibit  this  multiplicity  of  roots     and  therefore  do 
n 

not  indicate  higher  order  modes  of  operation  in  these  regions. 
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4.     Discussion. 

The  Ct!-j8  curves  for  the  finite  magnetic  fieid  case  show  several 
promising  facets.  The  backward  wave  mode  is  an  interesti  xg  pheno- 
menon, and  indicates  the  possibility  of  a  backward  wave  oscJ  la  tor. 

Fig.   (7)  shows  a  detailed  set  of  (jC  -  fi  relationships  for  the  case 
of  k  a  cot  0    =3,0,  b/a  =  0o8,  and  k  a  cot  $    =5.0.    Fig.    (8)  shows 
plots  of  the  left  hand  and  right  hand  sides  of  Eg.     3    3  5    vs    ka  cot  ij) 
for  pa  =  4.0.     The  perturbations  of  the  helix  wave  as  shown  in  Fig0   (7) 
indicate  that  considerable  coupling  exists  between  the  helix  and 
plasma  for  this  finite  magnetic  field  case.    This  Is  encouraging,  and 
indicates  that  one  might  reasonably  expect  jo  couple  the  input  and 
output  energy  directly  to  and  from  the  plasma  rather  thar  coupling  to 
the  electron  beam  and  thence  to  the  plasma  as  has  been  cone  by  previous 
workers,  and  this  direct  coupling  to  the  plasma  would  rhereby  allow 
excitation  of  the  backward  waves  as  previously  discussed. 

The  various  modes  of  operation  depicted  in  Fig.  e  plasma 

modes,  that  is,  they  are  not.  space  harmonic  modes  oi  the  helix, 
inasmuch  as  the  sheath  helix  approximation  does  no*      c  ude  the  space 
harmonics  which  exist  with  the  actual  helix.    Therefore,  in  a  solution 
which  utilized  a  more  exact  representation  for  the  hehx  one  might  ex- 
pect to  see   in  addition  a  coupling  between  the  plasma  and  higher  order 
helix  modes.     In  the  limiting  case  as  the  sheath  heiix  goes  over  to  a 
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conducting  cylinder,  these  modes  shown  would  be  the  various  radial 

modes  found  by  Trivelpiece,   (6)  and  denoted  as  the  p        modes  where 

ov 

p        is  the  l>th  zero  of  the  zero  order  Bessel  Function  of  the  first  kind. 
From  Fig.   (8),  we  see  that  the  solution  at  point  a  represents  the  per- 
turbed helix  wave,  the  solutions  to  the  right  of  point  a,  the  modes  in 
the  backward  wave  region,  and  those  solutions  to  the  left  of   a    the 
modes  in  the  forward  wave  region.    For  example  the  solution  at  point  b 
represents  the  first  plasma  wave,  and  this  point  corresponds  to  the 
point  b  on  Fig.   (7).    The  solution  at  point  e  represents  the  first  order 
solution  of  the  backward  wave  region  and  corresponds  to  the  point    e 
on  Fig.    (7).     The  other  points  correspond  as  they  are  lettered*     Thus 
we  may  see  the  mathematical  generation  of  these  higher  order  modes, 
and  may  visualize  them  as  representing  higher  order  radial  modes.     There 
are  of  course  an  infinite  set  of  these  modes,  only  a  few  of  which  have 
been  shown. 

The  slight  backward  wave  region  of  the  perturbed  helix  mode 
indicates  this  device  might  have  some  promise  as  a  backward  wave 
oscillator.     It  is  noted,  however,  that  the  bandwidth  of  any  such 
device  would  be  relatively  narrow,    The  forward  wave  characteristics 
indicate  that  a  forward  wave  amplifier  might  be  designed  to  be  quite 
broadband,  however. 

It  should  be  noted  that  the  higher  order  radial  modes  are  probably 
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not  particularly  useful  for  interaction  inasmuch  as  it  would  be  difficult 
if  not  impossible  to  couple  to  them  due  to  the  reversal  of  field  poten- 
tial in  the  radial  cross  section  which  such  a  mode  would  represent. 

Fig.    (9)  thru  Fig.   (16)  show  the  effect  of  variations  in  w    and  to 

P  c 

as  well  as  the  ratio  of  b/a.    These  curves  are  included  to  provide 

some  insight  into  the  effects  of  a  variation  of  these  parameters  .    First 

it  is  to  be  noted  that  for  oc    <  ic    the  lowest  order  plasma  wave   cuts 

pc 

off  at  05    ,  and  for  (£     <  co     the  plasma  wave  appears  to  cut    off  some- 
pep 

where  between  co    and  cc    .    The  backward  waves  and  helix  wave 
P  c 

approach  as  a  limit  ce    when  ui    >o>    ,  and  approach  as  a  limit  oe 

c  c        p  p 


when  Ul>    >  w   .    Further,  the  backward  waves  cut  off  at  thex  a;  s+  <x>  2. 
p         c  N    p  c 

The  above  observations  are  similar  to  those  found  by  Trivelpiece  (6) 

in  the  plasma  -  cylinder  problem.     The  larger  b/a,  the  more  completely 

the  plasma  fills  the  helix .     It  is  seen  that  this  leads  to  more  bandwidth 

in  the  backward  wave  region.    Also,  the  knee  of  the  curve  is  less 

sharp  as  the  helix  wave  goes  over  into  a  backward  wave  mode.    This 

one  might  expect  since  a  larger  b/a  would  lead  one  to  expect  tighter 

coupling  between  the  helix  and  the  plasma. 

Fig.    (17)  shows  a  plot  of  the  helix  and  plasma  waves  for  the  zero 

and  infinite  magnetic  field  case.    The  infinite  magnetic  field  case 

shows  a  continuation  of  the  effects  noted  previously.    There  appears  to 

be  a  smooth  transition  to  the  zero  magnetic  field  case  as  well.    As  cc 

c 
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is  decreased,  the  backward  wave  region,  essentially  between  a;    and 
\Jo;  a  +  05  2    becomes  smaller  and  smaller,  until  in  the  limiting  case 
there  is  no  such  region.    The  helix  wave  becomes  more  and  more  nearly- 
like  the  uncoupled  helix  wave  as  co     is  decreased,  and  finally  for  cc  =0 

c  c 

we  see  that  the  helix  wave  shows  no  backward  wave  region,  and  is 

only  slightly  perturbed  by  the  plasma  as  may  be  seen  by  comparison 

with  the  free  helix  wave  also  plotted  on  the  same  figure.    The  plasma 

wave,  which  at  values  of  uc    >  a;    appears  to  cut  off  at  oc    ,  cuts  off  at 

c         p  P 

lower  and  lower  frequencies  as  co     is  decreased  below  cc    and  finally 

c  p 

approaches  the  limit  as  shown  for  the  zero  magnetic  field  or  a;    =0 
case. 


22 


i 

i — [  j {-  4-  4- 

■  +  -fTffl 

ffWltff 

ff""  nnTnT"""" Tr :::  i :="~~ m rf  ^ 
ff  — rtttttii t~"  ----'*----  -------- 

r     1 III    III   ""  Tf T            ° 

f ~ ~  1  i II l"l  1  \ \  ~ ~    — 14- - ■               • 

r==l        M'^ifrr"           n 
n 

fc:  1              :±t:t+  :::             -a- 

1     1  [  1 1 1 1 1 1         i                     •«-• 

"  1— M ° 

""I                   ::                     ° 

-■-1  IIIIIIH '                    10 

a 

j-$|J 

o          '4  TO? 

CM 

II 

">        -a-         4-ffljr::::: 

ii        °        M|-- 

s.         o     -- 44114- 

•*       II    11   R  r  ' 

r  — 

l 

1 

| ; : ;   1 1 1 1 1 1 
^:fFj[r::{ 

^jjj---:]:llj_ 

|  j 1 1 1 1   1  j |f 

WW 

-LLi — 

pj-  mil  in  i  ii.i  i|i--ji-i — i — i — i — h-f  - 

IBhM1  ill 

ffl  ill       i  1  III  I  1  i  m# 

1  fV;  — 

■J — itHi  "    *■ 

11=     J 

r  III  {[I               n  t  fUffHt    * 

— — n — r i  rvH  ** 

1 — Lv^j 

D     :±: 
"*      ■+'■ 

m     tt: 

1 

J 

■ 

H 

B  1      jli  I--J- 

-  --     ll      - 'y-Jt-t- 

|  |  f[f^ H -i-j-l  -Li-4\  --- 

— ^ jii _ 

^4 

4 
( 
1 

a-     II 
»•      I1 1 

3     TTTTT4 

ffTrl 

TmTTTTT 

■■           > 

"--  ||   N4frT"f  """t""^~"  ++t^  J     r  ~*~~ 

^\         1 

::*?'«g;'ji|SI:::::: 

7T+4 

1 

2  Iee| 

Lfs  ;;;    ; 

.   ::. 

±i±±: 

::i::^H4|f:::::i-   ffl 

'      -■■.    '•}'■      ;     P\ 

iiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiiiiiiii 

j 

32 


5.    Conclusions. 

In  conclusion  it  appears  from  the  results  of  the  analysis  presented 
here  that  a  plasma  microwave  amplifier  device  in  which  the  R„  F.  Energy 
is  coupled  to  the  plasma  directly  offers  some  real  possibilities  of 
success.     It  appears  that  such  a  device  would  have  several  advantages 
over  conventional  type  traveling  wave  devices.     Principa  ly  these  are 
that  one  might  reasonable  expect  high  powe-  operation  in  a  small  device, 
and  that  one  might  be  able  to  alleviate  the  serious  mechanical  and 
cooling  problems  inherent  in  higher  frequency  miciowave  devices  of  the 
type  presently  in  use. 

However,  it  should  be  emphasized  that  the  analysis  presented 
contains  three  basic  approximations  which  limit  the  validity  of  the 
theory.     First,  the  quasi  -  static  approximation  is  made,  which  should 
be  good  for  the  low  values  of  beam  voltage  used  in  typical  tubes. 
Second,  the  sheath  helix  model  is  used,  which  is  believed  to  be  a 
reasonable  model  within  the  limitations  previously  mentioned.    Finally, 
the  assumption  of  a  free  space  TE  mode  within  the  helix  is  made.     This 
assumption  would  appear  to  be  justified  since  the  solution  obtained  in 
the  case  of  zero  magnetic  field  agrees  with  Paik's  analysis  (12)  which 
does  not  make  this  assumption. 

The  dispersion  relationships  of  Fig.     9    thru  Fig.    (16)  do  not  show 
any  anomalous  discontinuities  as  were  found  by  B.   M  Bulgakov  et  al. 
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(13). 

The  data  contained  in  this  paper  is  not  sufficient  in  itself  to  draw 
firm  conclusions  regarding  the  operation  of  a  plasma  microwave  ampli- 
fier device,  however,    The  0J-]8  diagram  does  not  provide  any  information 
regarding  the  magnitudes  of  fields  involved ,  and  thus  no  insight  into 
the  interaction  impedances  involved .    A  theoretical  study  of  interaction 
impedance  should  be  an  interesting  and  worthwhile  project  and  is 
suggested  as  the  subject  for  some  future  investigation. 

6.     Proposed  Experimental  Work, 

After  completing  all  connections  for  the  experimental  tube,  it  is 
proposed  that  the  following  experimental  work  be  accomplished i 

1.  Make  cold  tests  of  impedance  match  to  the  helix,  and  deter- 
mine the  losses  incurred  in  the  helix  couplers  over  a  frequency  range 
of  at  least  1.0  to  2,5  Gc. 

2.  Determine  the  CO  -  /8  relationships  as  indicated  from  the  theore- 
tical work. 

4.  Carry  out  electron  beam  interaction  experiments. 

5.  Investigate  the  backward  wave  regions  shown  theoretically, 
this  being  of  special  interest.     These  experiments  should  be  designed 
to  determine  the  range  and  magnitude  of  self  oscillation  available  as 
well  as  bandwidth  optimization  by  variation  of  field  strength  and 
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plasma  density. 

In  the  experiments  involving  measurements  of  the  R,  F.  interaction., 
measurements  should  also  be  made  to  determine  the  impedance  jevels, 
gain,  or  loss  in  the  helix  regions,  and  gain  or  loss  in  the  drift  region. 
Further,  these  experiments  should  be  carried  on  over  a  wide  range  of 
magnetic  fields,  tube  temperatures,  and  frequencies, 

Care  should  be  exercised  to  assure  the  cooling  provided  by  the 
water  jacket  is  adequate  to  maintain  an  effective  temperature  along 
the  region  so  cooled.    Fig.   (18)  shows  a  photograph  of  the  laboratory 
equipment  designed  to  work  with  the  proposed  tube.    The  magnet,  tube, 
and  heat  sinks  for  the  filament  rectifiers  are  water  cooled,    fig.   (19) 
shows  a  block  diagram  of  the  proposed  measurement  setup 

The  magnet  is  designed  to  provide  an  adjustable  uniform  field 
from  0  to  2000  gauss.    The  plasma  generation  filament  supplies  are 
designed  to  provide  from  two  to  six  volts,  and  up  to  a  maximum  of 
120  amps  each.    The  rectifiers  are  rated  at  250  amps,  ,  the  transformer 
being  the  limiting  element. 
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APPENDIX  I 
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..JOB  JONES  GRAPH  CM  8   7 
PROGRAM  jn^ESS 
PRINT  1 

1  FORMAT  (12H  TUBE  DESIGN   ) 

16  FORMATM8H  IMAGINARY  GAMMA     ) 
15    FORMAT  (2X,6(2X,E15.8) ) 

5    FORMAT  (5x,6(2X,E15.0n 

17  FORMAT(  18H  X  ,  Y  ,  HZ  ,  HZ  Y,  F,-HBB    ) 
7     FORMAT!  5X  ,  3  (5X,C  15.8)  ) 

PRINT     17 

DIMENSION    XI (900) ,Y1(903 ) 

DIMENSION    T( 1000) 
C  Y=    KACO+    PSI 

C  X=    PA 

C  C0N1     =KPAC0T(PS1)      . 

C  C0N2=B/A 

C  F=    -J(OMEGA) (EPS010N    ZERO)ZP(YS) 

CON3=5.0 

C0N2=0.8 

CON1=3.0 

PRINT  7.C0N1 ,C0N2,C0N3 

11=50 

LL=200 
43         CONTINUE 

DO   40    M=LL,1 500C, II 

L  =  M 

KA=-1 

3  DO    50    K=5,375, 1 
53         X=FLOATF(K)/25.0 

N  =  L 

AA=C0N2*X 

AO=BESI (0.0,X) 

A1=BESI ( 1.0, X) 

CALL    PESS(0,X,2,B0,T) 

CALL    HESS( 1, X,2,B1 ,T) 

CO  =  BESI  (CO,  AA) 

C1=RESI( 1.0,AA) 

CALL    RESS(0, AA,2,D0,T) 

CALL    PESS( 1, AA,2,D1,T) 

4  DO    60    J=N,     15000,     1 
Y=FLOATF( J)/2C00.0 
L  =  J 

A=. 0000001 

IF(SENSE    SWITCH    2)     102,103 

102  Y  =  2.«Y1 (KK-1  )-Yl (KK-21-A 

103  AB=Y/X 
Q=(C0N1/Y) 
RA=Y**2.0 
RB=CON3»*2.0 
RC=RA-RB 

51  RD=(C0N1 )»*2.0/RC-1.0 

RR=1  .0-(CONl  )**2.0/RC 

GAP    =    ABSF(( UQ)*«2.0)-1.0)/RD)*«0.5 

GAB    =    C0N2*X*    GAP 
-      XNUM1    =(AB)*«    2.0 

DEN1=A0/A1     +B0/B1    -     ( B 1/ BO ) »XNUM 1 

F    =XNUM1/DEN1 

IF    (SENSE    SWITCH    1)     101,120 
120       IF(Q-l.O)     12,13,13 

12       IF(RD)     10,11,11 
13  IF    (RD)     1  1  ,99,10 

99         PRINT    2 

2  FORMAT     (16H    GAMMA     INFINITE  ) 
11          CALL    BES( 1,GAB,0,BESJ1,T ) 

CALL    BES(0,GAB,0,BESJO,T) 
E1=BESJ1 
EO    =BESJO 
BB=RR«GAP*E1/E0 
91  CC=C1*(B0+F*(B1 ) ) 

DD=    Dl*     (AO-    F*( A1 )) 
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EE=-CC»(B0+F*(B1  )  ) 

TF=    DOMAO-FM  Al ) ) 

XNUM2    =    CC+DD 

DEN2    =    EE    +FF 

Z=XNUK2/       QEN2 

ZY=Z/BB 

ALFA-ACSF (ZY-1. ) 

EPS=0.1 

IF     (ALFA-EPS)    20,20,59 

20  PRINT    5,X, Y, Z,ZY,F,BB 
GO    TO   95 

10         HAP    =    -GAP 

FO    =    PESKO.CGABJ 

Fl=     BFSH  1.0,  GAB) 
93         HCC=C1»(BC+F»(B1 ) ) 

HDO=    Dl*( AO-F*( 41 )  ) 

HEE    =    -C0*(B0+F*(B1 ) ) 

HFF=    00*( AO-F* (A1  )  ) 

HB8=RR«HAP*F1/FC 

DHN=HEE+HFF 

HNM=HCC+HDD 

HZ=HNM/DHN 

HZY=HZ/HPB 

HLFA=    ABSF     (HZY-1.0) 
•     EPS=0.1 

IF(HLFA-EPS)21,21,59 
59  IFCSENSE    SWITCH    2)     105,60 

105  1F(SENSE  SWITCH  3)106,107 

106  A=A+. 000005 
GO  TO  102 

107  A=A-. 000005 
GO  TO  102 

60   CONTINUE 

21  PRINT  15,X,Y,HZ,HZY,F,HBB 
PRINT  16 

95    KK-K-l* 

XI (KK)=X 

Yl (KK)=Y 

NUMPTS=KK 
50  CONTINUE 

101       PRINT    5,X, Y,Z,ZY,F,BB 

PRINT    5,HZ,HZY,r  ,HBe.Q,R() 

CALL    GRAPH    ( NUMPTS ,X 1 , Yl , 8 ) 

IF(SENSE    SWITCH    3)     12.8, U  1 
128       11=100 

LL=M+I  I 

GO  TO  43 
41  CONTINUE 
40    CONTINUE 

STOP 

END 

SUBROUTINE  BESS  ( NO, X, KOOE, RESULT, T > 

DIMENSION  T(  1000  ) 

PI=3. 1415926536 

IF  (KOOE)  41,1,2 

1  CALL  PES  (NO, X,0, RESULT, T) 
RETURN 

2  IF  (KODE-1  )  41  o,4 

3  CALL  BES  ( NO, X , 1 , RESULT, T ) 
RETURN 

4  IF  (KODE-2)  41.5.6 

5  IF(N0-1 )   11,11,12 

12  PRINT  13 

13  FORMAT  (34HlORDtR  ERROR  FOR  K  AND  Y  FUNCTIONS) 

14  STOP  14 

11  IFCX-1.1  15,15,16 
16  T( 1 )=1./X 

DO  240  1=2,  12 
240  TCI )  =  T( 1-1  )»T<1) 
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IF(NO)   21 

20  XG0=1. 253.3 
1 T(3)+.134L 
2-. 52472773 
3+.218M5180 

RESULT  =  E 
RETURN 

21  IF  (NC-1) 

22  XG1  =  1.253 
1+.  128014266 
2-.459M3421 
3+. 50502385 
4-. 01082417 

RESULT  =  E 
RETURN 
15  CALL  RES  ( 
XIO=T( 1 ) 
T(2)=X*X/4 
DO  310  I=U 
310  T(  I  )=T(  1-2 
SUM 


2  C  2  1 
14  13"  73-.  156664  1  8  1  6 *  T  {  1  J+.C881  1 
59  622lJ*T(4  )-.2  299850328*T(5)+. 
31 *T( 7 ) +.5  5755  68  367*1 (8)-. 4262 
96*T< 10)-.0  66B097672*T( 11 J+.00 
XPF(-',  .*X)/SG*TF(X)«XG0 


12,22,1 
3  1  4 1 37  3 
36*T(3)- 
17*T(6)+ 
76*T( 9)- 
75«T( 12) 
XPF(-1 .* 


+.U699927013*T( 1  J-.146 
.1736i431637*T(U)  +  .28U7 
.6  28 33  80  68 1*T( 7)-. 66 32 
.2581303765*U  10)  +  . 078 


12  7P2*T(2)-.09  131095U6« 
3792U09730M (6) 
632912»T (9) 
91B93330*T{ 12) 


8582957*T(2) 
6181 49*T(5) 
295U30»T(8) 
8000113*T{ 11 ) 


X)/SQ^TF(X)«XG1 
1, X, 1 , XI 1,T) 


-.5 

+  .0 

2        +.0 

ANSKO  =  SU 

IF(NO)   31 

30  RESULT  =  A 
RETURN 

31  IF (NC-1) 

32  RESULT  =  ( 
RETURN 

6  IF  (KODE-3 
Ul  PRINT  42 

42  FORMAT 

43  STOP  43 
140  IF  (NO-1) 

50  T(l)=X/3. 
DO  UOO  1=2 

UOO  T(  I)=T(  1-1 
IF  (3.-X) 

51  SUMY  =   .3 

1  +.2 

2  -.0 
CALL  RESH 
XJO=T( 1 ) 
Y0=2./PI*L 
IF  (NO)   5 

54  RESULT  =Y0 
RETURN 

53  IF(N0  -1 ) 

55  RESULT=(XJ 
RETURN 

52  IF  (  4.-X) 
60  T(  1  )=X/3. 

DO  402  1=2 
402  T(  I  )=T( 1-1 
SUMFO  =.79 
-.00 


,12,2 
)*T(2) 
7721566 
34  88590* 
0000740* 
M  -  LCGF 
,30,31 
NSKO 

12,32, 12 
l./X  -AN 

)   U 1  ,  14  0 

(1  1H1K0D 

5  0,50,1 

,12 

)*T(  1) 

52,51,51 

6746691 

5300117* 

0024846* 

,X,0,XJ1 

0GF(X/2. 
3,54,5  3 


+.U2278420 
T(6)+.  00262698 
T(  12) 
(X/2.)  *  XIO 


SKO*XI  1  )/XIO 
,41 

E    ERROR) 
2 


+  .60559366  * 
T(6)-  .04  261214  * 
T(  12) 
,T) 

)*XJO+SUMY 


*T(2)   +  .23069756  *  T(4 ) 
*T(8)   +  .00010750  *  T( 10) 


T(2)  -  .74350384  *  T(4) 
T(8)  +  .00427916  »  T( 10) 


12,55,12 
1*YC-2./(PI*X) )/XJO 

59,60,60 


1 
2 

1 

2 


70 

71 
555 


59 


+  .00 
SUMS   =.78 
-.00 
-.00 
YO  =  SUMFO 
IF  (NO) 
RESULT  =Y0 
RETURN 
IF  (NO-1) 
CALL  PES( 1 
XJO=T( 1 ) 
GO  TO  55 
T(l  )=4./X 


,12 

)«T(  1  ) 

78845ft 

009512/T 

C14476/T 

539816 

262573/T 

013558/T 

/  SQRTF( 

71 ,70,71 


12,555,1 
,X,0,XJ1 


-.D0CC0077  #3./X 
(3J+.D0137237  /T(4) 
(6) 

+.04166397  *3./X 
(3)+. 00054125  /T(4) 
(6) 
X)  *  S INF  (X-SUMS) 


2 
,T) 


-.00552740/1(2) 
-.00072805/  T(5) 

+.00003954/  T(2) 
+.00029333/  T(5) 
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DO  401  1=2,12 
UOl  T(I )=T( 1-1 )*T( 1 ) 
IF  (NO)   12,73, 12 

73  PO  =  .3989U2279:5      -.0017530620  *T(2)   +.0001734300  *r{4) 

1     -.000048761;>*T(6)  +  . 30001  73565  *T(8)   -.OO00C370U3  »T(  10) 

QO  =   -.012U669UU1      + .000456U32U«T(2)   -.0000869791  *  T(U> 
1       ♦.000O5U2468»T{6)-«0nooiM?078«T(fl)   ♦.00000323  12  «  TMO) 

YO  =  2./SQRTF(X)*(P0*SIMF(X-PI/4. ) +  U . /X*GO*COSF ( X- PI /4. ) ) 

RESULT=YO 

RETURN 
72  IF  (NO-1)    12,7U, 12 

74  PI  =     .3989^22819     +  .00292  1  8256*  T  (  2  )    -.0002232030  *T(U) 

1       +.0000580759«T(6)-.C000200920*T(8)   +.0000042414  *  T(10) 

Ql  =    .0374008364      -  .0006390UC0»T(  2  )   +.0001064741  »  T  { 4  ) 
1       -.0000398708*T(6)f . 00001 62200*T ( 8 )   -.000003659U  *  T(10) 

Yl  =  2./SQRTF(X)«(Pl»SI\iF(X-3.«PI/4.  ) +4. /X*Q 1 *COSF ( X-3. *P  1/4.  )  ) 

RESULT=Y1 

RETURN 

END 

SUBROUTINE  B  ES  ( .MO,  X,  KODE  ,  RESULT  ,  T  ) 

DIMENSION  T( 100) 

KLAM=1 

K0=N0+1 
1009  IF(X)  81,101,82 

E  RR= • CO  1 

81  IF'  (  X+ERR) 108, 101 •  101 

82  IF(X-ERR) 101,101 , 108 

101  IF(NO) 103, 102,  105 

102  T(K0)=1.0 
RESULT  =1.0 
RETURN 

103  IF(KO) 106, 104, 105 

104  RESULT=9.999999999E200 
RETURN 

105  RESULT=0.0 
RETURN 

106  PRINT  107 

107  F0RMAT(55U1NEGATIVE  ORDER  NOT  ACCEPTED  IN  BESSEL  FUNCTION!  ROUTINE) 
1071  STOP  1071 

108  IF(NO) 106, 1081,1081 
1081  IF (KODE) 1 1  ,12,11 

11    KLAM=KLAM+1 

12  J0=2*XABSF(XFIXF (X) ) 
MO*  NO 

IF(MO-JO)2,21,21 
2     MO  =  JO 

21  M0=M0+11 

22  T(M0)=0. 
LUB=MC-1 
T(LUB)=1 .0E-300 

GO  TO  (23,51  )  ,KLAM 

23  F=2*LUB 

231  ENA(-3)RAD(MC) 
LIL2(M0) ,ENA(T). 
INA( 1 )SAU(233P) 
INA(  1  )SAU(233A) 
INA(  1  )SAL(233A) 

232  LDA(F)FSB(2. )  ' 
STA(F)FDV(X) 

233A  FMU2(N)FSB2(N) 
233B  STA2(N) IJP2(232) 

SUM=T( 1 ) 

GO  TO  (39,69),KLAM 
39  M3=M0+3 

DO  4  J=3,M3.2 

4  SUM  =  SUM+2.*T(  J) 
F=l ./SUM 

DO  5  J=1,K0 

5  T(J)=T(J)*F 
RESULT=T(KO) 
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51 
69 
7 

8 


301 
302 

202 

203 
205 

20U 
206 

208 

209 

211 
213 

212 
216 

214 

215 


M 

6 


RETU 
K0N  = 
LDA( 
GO  T 
LAC( 
M3«M 
DO  7 

S 
F  =  l. 
DO  3 

T 
RESU 
END 
SUBR 
DIME 
.6(1) 
6(2 
B(  3 
B(4 
B(5 
B(6 
B(7 
B(8) 
B(9) 
B(  10 
B(  1  1 
B(12 
B(  13 
Z  = 
I  = 
z  = 
I  = 

IF  ( 
Z  = 

END 


RN 
777 
KCN 
0  2 
KCN 
0+3 
J» 
UM  = 

/su 

J  = 

(J) 

LT  = 


7777776777777B 
),RAD(235A) 

) , RAD (233  A) 

2,M3 

SUM+2.*T( J) 
M*EXPF(X) 
1  ,K0 
=T( J)*F 
T(KO) 


CU 
NS 

)  ■ 

)  = 

)  = 

)  = 

)  = 

)  = 


) 

) 

) 

) 

B( 

12 

Z* 

I- 

1  + 

(Z 


TINE  R 

ION  8( 

-.U22 

-.233 

.1910 

-.02U 

-.017 

.0080 

-.000 

-.000 

.0001 

=  -.00 

=  -.00 

=  .000 

=  -.00 

13) 


GAM 

13) 

7845 

C937 

91  1C 

5524 

6452 

2  32  7 

804S 

3608 

4562 

0017 

C002 

con 

0000 


(X,Z) 

35092 

36565 

1162 

908°7 

42118 

81  13 

41  335 

51491 

4324 

527917 

625721 

23554 

181220 


X  +  8(  I) 
1 

1  ) 
+  1 


302,302,301 
)»X*(X+1.  ) 


SUBROUTINE  RCPGAMR(X,Z) 
IF  (X-l.)  208,202,203 
Z  =  1. 
RETURN 

Y  =  1. 

X  =  X  -  1  . 

Y  =  Y«X 

IF  (X-l.)  206,204,205 

Z  =  l./Y 

RETURN 

CALL  RGAM  (X,Z) 

Z  =  Z/Y 

RETURN 

INT  =  X 

Q  =  INT 

IF(X-Q)  211,209,211 

Z  =  0.0 

RETURN 

IF  (X+l.)  212,209,213 

CALL  RGAM(X.Z) 

RETURN 

Y  =  X 

X  =  X  +  1  . 

IF     (X+l.)     214,209,215 

Y  =    Y*X 
GO    TO    216 
CALL    RGAM(X,Z) 
Z    =    Y*Z 

END 

FUNCTION    BESI(FN,Z) 

IF(Z)    4,5,5 

PRINT    6  . 

FORMATt     U2H    BESI     NOT     PROGRAMMED    FOR    NEGATIVE    ARGUMEN1 

PRINT    7 
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7 
5 

9 

8 


13 
12 


100 
101 


102 
105 


103 


104 
200 
201 
202 

203 

205 


204 
300 


FORM 
Z  =  AB 
•  I  F  (  M 
FN  =  A 
FK=1 
FNl- 
FN2» 
CALL 
CALL 
AI  =  ( 
A  =  ZA 
FK  =  F 
AI  =  A 
FN2  = 
A  =  A  + 
IF(A 
BESI 
END 
FUNC 
FN  =  A 
IF(M 
FA=- 
DESK 
RETU 
N  =  FN 
A={- 
HL  =  L 
BIBA 
BIA  = 
B1B= 
IF(N 
b  =  0. 
GO  T 
DO  1 
FKK  = 
BIBA 
BIB  = 
BI  =  B 
B  =  BI 
FK=1 
FNK  = 
BIA  = 
BIB- 
BI=B 
FK=F 
TEST 
B  =  B  + 
IF(A 
B=(- 
IF(N 
C  =  0. 
GO  T 
C=l. 
GO  T 
NA=N 
NN=1 
DO  2 
NN  =  N 
C=NN 
CB=C 
KD=1 
FKD  = 
CB=( 
KD=K 
C=C  + 
IF(K 
C  =  C* 
BESK 
END 
END 


AT(  U5H  BESI  COMPUTED  FOR  ABSOLUTE  VALUE  OF  ARGUMCNT     ) 

SF(Z) 

CDF(FN,1  .)  )    8,9,8 

BSF(FN) 

FN+1. 
FN*ii 

RCPGAMR(FN1 ,  ZA> 

RCPGAMR(FN2,ZB) 
1/2. )**2»ZB 
+A  I 
K+l  . 

I*( Z/2. )**2/ ( (FK+FN) *FK) 
FN2+1. 
A  I 

BSF( A  I  )-ABSF(A)*lO.*«(-14) )  2,3,3 
=(Z/2. )**FN*A 

TION  BESK(FN,Z) 

BSF(FN) 

0DF(FN,1.))  13,12,13 

FN 

=  1. 57079 6327 »(3ES  I  (F A,Z)-BESI ( FN , Z ) ) / ( S I NF ( 3. 1 41  592 65 4* FN  )  ) 

RN 

!.)**( N+1)»BESI(FN,Z )*LOGF( 1 .78 1 072418*Z/2 . ) 

/2. 

=  1. 

0. 

0. 

)  100,100,101 

0  105 

02  KK=1,N 
KK 

=BIBA/FKK 
BIB+1./FKK 
IBA 
BA*BIB 

FN  +  FK 

BIA+l./FK 

BIB+l./FNK 

I«HZ**2/(FK*FNK) 

K+l  . 

=BI*(BIA+BIB) 

TEST 

BSF(TEST)-AHSF(B  )  » ID . »» ( - 1 4 ) )  104, 103, 103 

HZ  )**N«B 

-1)  200,201,202 

0  300 

0  204 
-1 

03  J=1,NA 
N*J 


KD 

-1. )«*KD*CB*HZ«*2/U  FN-FKD)«FKD) 

D+l 

CB 

D-N)  205,204,205 

HZ**(-N) 

=A+(B+C)/2. 
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0 

02 
02 
0 


02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 

0 

02 

02 

0 
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OMEGA  BETA  PLOT ,C 1  =  3 . , C2  =  0.  b  C3=2.0 
COLIN  M.  JO.NCS 


OMEGA  BETA  PLOT ,C 1»4. ,C2«0.  5  CS«2.0 
COLIN  H.  JONES 


OMEGA  BETA  PLpT ,C 1 =3 . ,C2  =  0.  b  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3 . ,C2  =  0.  5  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3 . ,C2  =  0.  b  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3 . , C2  =  0.  b  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT »C 1 =3 . ,C2=0.  b  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1=3 . ,C2=0. 5  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3. ,C2  =  0.  5  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1=3. , C2  =  0. 5  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1  =  3. ,C2  =  0.  5  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1  =  3. ,C2=0.  b  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1  =  3. ,C2  =  0.  5  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1  =  3 . ,C2  =  0.  b  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3. ,C2=0.  b  C3  =  2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT ,C 1 =3 . ,C2  =  0.  b  C3=2.0 
COLIN  M.  JONES 


OMEGA  BETA  PLOT , C 1 =3 . ,C2  =  0.  b  C3=2.0 
COLIN  M.  JONES 


10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
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